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Abstract. We report on the measurement of a Larmor frequency shift proportional to the electric-field strength for 
^®®Hg atoms contained in a volume permeated with aligned magnetic and electric fields. This shift arises from the in¬ 
terplay between the inevitable magnetic field gradients and the motional magnetic field. The proportionality to electric- 
field strength makes it apparently similar to an electric dipole moment (EDM) signal, although unlike an EDM this 
effect is P- and T-conserving. We have used a neutron magnetic resonance EDM spectrometer, featuring a mercury 
co-magnetometer and an array of external cesium magnetometers, to measure the shift as a function of the applied 
magnetic field gradient. Our results are in good agreement with theoretical expectations. 

PACS. XX.XX.XX No PACS code given 


1 Introduction 

Recent investigations characterizing frequency shifts for spins 
contained in vessels permeated with magnetic and electric fields 
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B, E have been motivated principally by the search for elec¬ 
tric dipole moments (EDMs) of simple non-degenerate sys¬ 
tems (neutron, atoms, molecules) and the potential discovery 
of new sources of CP violation ||T|. Such experiments look for 
shifts, proportional to an applied electric field, of the Larmor 
precession frequency of stored particles. Any additional such 
shift is therefore a potential source of systematic errors. Among 
the few magnetic-field related spurious shifts, one is of par¬ 
ticular concern: due to the motional magnetic field E x v/c^, 
a shift arises that is proportional to the electric-field strength 
and therefore mimics an EDM signal. Interestingly enough, 
E X v/c^ effects were already the main limiting factor for the 
early neutron beam experiments E). Then, with the advent of 
the storable ultra-cold neutrons (UCN), it was erroneously as¬ 
sumed for many years that this false EDM signal would vanish, 
based on the argument that the velocity of trapped particles av- 
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erages to zero. The first correct and comprehensive calculation 
of this effect was given in Ref. 0, in the context of an EDM 
experiment with stored particles. For completeness, it should 
be mentioned that Stark interference effects, such as the one 
reported for in Ref. |]4l, are also known to produce false 

EDM signals for atoms. The effect discussed in the present ar¬ 
ticle is of a different nature, and to make the distinction we will 
refer to it as the motional false EDM. 

Our collaboration is conducting a program to search for 
the neutron EDM Q, using the new ultracold neutron (UCN) 
source 0 at the Paul Scherrer Institute (PSI). We are currently 
working with an upgraded version of the spectrometer Q that 
was used to establish the best nEDM limit, 

\dn\ < 2.9 X 10-2®ecm(90%C.L.), 

at the Laue Langevin Institute (ILL) jS). One distinct feature 
of this device is a mercury co-magnetometer 0 using a spin- 
polarized vapor of ^®®Hg atoms that precess in the same vol¬ 
ume as the neutrons. The nEDM analysis is then based on the 
ratio of the Larmor precession frequencies, R = fn/fag, which 
to first order is free of magnetic field fluctuations. However, 
both neutrons and mercury atoms are subject to a frequency 
shift that is proportional to the electric field, due to the unavoid¬ 
able presence of magnetic-field gradients. As will be shown, 
the motional false neutron EDM, is negligible, at least at 
the current level of sensitivity. In contrast, the mercury-induced 
false nEDM 

^false.Hg ^ _ 3 g (j) 

7Hg 

where ® is the motional mercury false EDM and 7n, 7Hg are 
the gyromagnetic ratios of the neutron and ^®®Hg respectively, 
is a major systematic effect that must be precisely controlled. 

One of the main improvements accomplished recently within 
the experiment is the installation of an array of cesium mag¬ 
netometers that surrounds the precession chamber. This new 
device has made it possible to measure the magnetic field dis¬ 
tribution, and thus to calculate the vertical gradient in the trap, 
which underlies the false EDM discussed here. 

In this article, we report on the first direct measurement 
of a motional false EDM signal for stored mercury atoms. A 
comparison to theoretical expectations is also presented. 

2 Theory of frequency shifts induced by 
magnetic fieid gradients: a brief reminder 

Particles with a magnetic moment exposed to a magnetic field. 
Bo = BqZ, precess at the Larmor frequency /l = ^ Bq/2tt 
where 7 is the gyromagnetic ratio. Because of experimentally 
unavoidable magnetic field gradients, the Larmor frequency of 
a particle moving through this field will be subject to a shift, 
known as the Ramsey-Bloch-Siegert (RBS) shift flO). If an 
electric field E (parallel or anti-parallel to Bq) is applied - as 
is the case in experiments searching for EDMs - the moving 
particle will experience an additional motional magnetic field 
B„ = E X v/c2. It is the interplay between this field and the 
magnetic field gradients that lies at the origin of a frequency 


shift proportional to the electric field strength, thus inducing a 
false EDM. 

As mentioned above, the first detailed calculation of such 
false EDMs for stored particles was given in Ref. a in the 
context of the RAL-Sussex-ILL neutron EDM experiment 0. 
The authors derived expressions for the two limiting cases: 
non adiabatic and adiabatic, corresponding to 27r ^ 1 and 

‘C 1 respectively, where t is the typical time parti¬ 
cles take to cross the trap. Both regimes are of interest, since 
is^Hg atoms fall into the first category whereas UCNs fall into 
the second. More general results, valid for a broad range of 
frequencies, were obtained only for cylindrical symmetry and 
specular reflections. The expressions of the frequency shifts for 
the two limiting regimes are : 

^2jj2 

^/l = 17 ,—(non adiabatic) (2) 

SZTT OZ 

(adiabatic), (3) 

47r Bq & OZ 

where 7 is the gyromagnetic ratio, D is the diameter of the 
trap, c is the velocity of light and Vxy is the particle velocity 
transverse to Bq. Note the absence of the gyromagnetic ratio in 
Eq. 0. Indeed, in the adiabatic case, the frequency shift can 
be interpreted as originating from a phase of purely geometric 
nature, or Berry’s phase mm, and is therefore independent 
of the coupling strength to the magnetic field. 

These results were then complemented and extended us¬ 
ing the general theory of relaxation (Redfield theory) 013II14I . 
and then by solving the Schrodinger equation directly ifBU . In 
Ref. m, an expression valid for arbitrary field distributions or 
trap shapes was obtained in the non-adiabatic limit: 

7 ^ 

'^/l = — k {xBx + yBy) E (non adiabatic), (4) 

27rc^ 

where the brackets refer to the average over the storage volume. 
For a cylindrical uniform gradient and a trap with cylindrical 
symmetry, Eq. 0 reduces to Eq. •Ill- 

Using the relationship between the frequency shift and the 
false EDM, 

where h is Planck’s constant, together with Eqs. 0 and 0, one 
can now readily calculate the magnitude of the false EDMs for 
the mercury and for the neutron (both direct and mercury in¬ 
duced). Given our experimental conditions (see section]^ and 
assuming a neutron velocity of 3 m/s, one obtains: 

Q D 

^faise ^ 0 ^ ^ 10-29 ecm/(pT/cm) (5) 

OZ 

f) n 

1.148 X 10-27 ecm/(pT/cm) (6) 

no 

rfaise.Hg ^ 0 ^ ^ q -27 ^ cm/(pT/cm). (7) 

OZ 
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Considering a typical value of 10 pT/cm for the vertical (z di¬ 
rection) gradient in our setup, we can conclude on the one hand 
that the direct false neutron EDM is negligible, at least at the 
current level of sensitivity. On the other hand, the mercury- 
induced false neutron EDM is a major systematic error that 
must be properly taken into account. 


3 Experimental apparatus 


The experimental study was performed with the nEDM spec¬ 
trometer installed at the PSIUCN source. This room-temperature 
apparatus uses the Ramsey method of separated oscillatory fields 
I MT^ to search for a shift, proportional to the strength of an ap¬ 
plied electric field, in the neutron Larmor precession frequency. 

Under normal operation, polarized UCNs are stored in a ~ 
20 liter chamber (internal diameter D = 47 cm, height//= 12 cm), 
composed of a hollow polystyrene cylinder (coated with deuter- 
ated polystyrene) 11711181 and two disk-shaped aluminum elec¬ 
trodes coated with diamond-like carbon (Eig. [T]). A cos0 coil 
produces a highly homogeneous magnetic field, Bq ~ 1 pT, 
in the vertical direction while the two electrodes - the top one 
being connected to a high voltage (HV) source and the bot¬ 
tom one to ground potential - generate a strong electric field 
(E « lOkV/cm), either parallel or anti-parallel to Bq. In addi¬ 
tion, a set of trim coils permits an optimization of the magnetic 
field uniformity at the 10 “^ level. 


Top CsM 

25.1cm 


Hg Lamp 


Bottom CsM 
Super bottom CsM 


A .a. 



' Upper electrode 
{HV) 


PMT 

• Lower electrode 
(Ground) 

. Hg polarization 
chamber 


Fig. 1. Schematic view of the precession chamber of the nEDM @ PSI 
experiment. 


The key to such experiments relies on the ability to con¬ 
trol the magnetic field both in terms of stability and homo¬ 
geneity. To this end, we use two highly sensitive and com¬ 
plementary atomic magnetometers based on mercury (^®®Hg) 
and cesium (^^^Cs) atoms, respectively. Mercury is used in a 
co-magnetometer mode: polarized mercury atoms precess in 
the same volume as the neutrons, hence probing approximately 
the same space- and time-averaged magnetic field. Cesium is 
used in a set of external magnetometers surrounding the stor¬ 
age chamber. The former is an ideal tool to correct for field 
drifts, while the latter gives access to the spatial field distribu¬ 
tion. 


3.1 The mercury co-magnetometer 

To date, is the only atomic element that has been used 

as a co-magnetometer for a neutron EDM experiment. Thanks 
to its nuclear polarization, it benefits from long wall collision 
relaxation times, and polarization lifetimes larger than 100 s 
can be achieved. Moreover, it is one of the rare elements in 
which nuclear spin polarization can be created and monitored 
by optical means. It is worth noting that the best absolute EDM 
limit comes from an experiment using ^®®Hg|| 22 |Q 

|(i(i®®Hg)| < 3.1 X 10-29 ecm (95%CL). 

In our experiment, a vapor of mercury atoms is spin-polarized 
by optical pumping in a polarization chamber located under¬ 
neath the precession chamber (Fig.[^. The operation of the co¬ 
magnetometer is synchronous with the nEDM measurement, 
and follows cycles about 300 s long. During neutron counting 
and filling, mercury atoms are continuously injected and opti¬ 
cally pumped in the polarization chamber. Once the precession 
chamber is filled with UCNs, we let the vapor diffuse into the 
precession chamber where, after the application of a 7r/2 pulse, 
the atoms freely precess around Bq at a frequency of about 
8 Hz. The interaction of the precessing atoms with a circu¬ 
larly polarized resonant probe beam produces a light-intensity 
modulation whose analysis yields the Larmor frequency of the 
atoms. 

One of the major drawbacks of the ^99 jjg co-magnetometer 
is its sensitivity to high voltage. As illustrated in Fig.|^ which 
displays the transverse polarization relaxation time T 2 versus 
the cycle number, sudden T 2 drops are systematically observed 
after each HV polarity reversal. The corresponding reduction of 
the signal amplitude directly affects the precision of the mag¬ 
netometer. Fortunately, optimal performance can be recovered 
via discharge cleaning in an oxygen atmosphere. On average, 
the precision of the mercury co-magnetometer is of the order of 
100 fT, equivalent to a magnetometric precision at the 0.1 ppm 
level per cycle. 


3.2 The array of cesium magnetometers 

An array of 16 cesium magnetometers (CsM) m allows mea¬ 
surement of the magnetic field distribution in the region of 
interest and, in particular, it gives us knowledge of the verti¬ 
cal gradient dB^Idz. Six HV-compatible (i.e. fully optically 
coupled) magnetometers were placed on top of the precession 
chamber, and ten standard ones below (Fig. [T| and [^. These 
laser-pumped magnetometers use a vapor of ^Cs atoms (gy- 
romagnetic ratio 7 = 27r x 3.5kHz/pT) and are operated in a 
phase-stabilized mode. They have a high statistical sensitivity 
(^ 100 fT for 40 s long measurements); however, they suffer 
from inaccurracies of their absolute field readings, with offsets 
that can be as high as 100 pT. They are therefore precise but 

* One may wonder why the effect discussed in the present article 
was not observed in that experiment. They actually use spectroscopy 
cells filled with 475 Torr of CO buffer gas acting as a UV quencher. 
Consequently, mercury atoms move in the diffusive regime where the 
motional false EDM essentially vanishes - in contrast to the ballistic 
regime of our mercury co-magnetometer. 
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Fig. 2. Transverse relaxation time T 2 of atoms (green points) 

together with the high voltage value (blue line) versus cycle number. 
Sudden drops of T 2 are observed after each polarity reversal. 


not accurate. Finally, it is important to note that these magne¬ 
tometers - like the mercury co-magnetometer - are scalar: they 
measure the magnitude of the magnetic field at the center of the 
bulb containing the cesium vapor. 



Fig. 3. Picture of the six HV-compatible Cs magnetometers installed 
on the top HV electrode in A1 enclosures. Optical fibers are also visi¬ 
ble. 


4 Measurement and data analysis 

A preliminary measurement with a limited number of CsM was 
performed in 2011, and led to a first result EOll . The present 
analysis is based on a dedicated data-taking period of 2 weeks’ 


duration in December 2013, where eight different gradient set¬ 
tings were explored: four with the magnetic field pointing up¬ 
wards (iJp), and four downwards (Bq). Two trim coils were 
used to set a vertical gradient in addition to the Bq field gener¬ 
ated by the main coil. For each field configuration, about 500 
cycles were recorded with a basic FIV polarity pattern (-1- — 
—h) and polarity changes every 20 cycles. The voltage was set 
to 120 kV, i.e. as high as possible to maximize the frequency 
shift while preserving a smooth operation (limited number of 
electrical breakdowns). 

As discussed above, the frequent polarity reversals induced 
a significant degradation of the mercury magnetometer’s sen¬ 
sitivity. Consequently, we decided to limit the free precession 
time to 40 s, a good compromise between sensitivity and the 
number of cycles. The mercury frequency was extracted using 
our standard “two windows” method ETll . It consists of fitting 
the signal phase at either end of the signal, using data in two 
15 s windows at the beginning and end of the time series. This 
method optimally takes into account possible frequency drifts 
during the precession time. During data taking, the mercury 
frequency uncertainty varied in the range 1-2 /iHz. 

Outputs from all 16 CsM were continuously recorded at 
a rate of 1 Hz, and a mean value of the magnetic field was 
calculated for time periods having an exact overlap with the 
mercury precession. We further made the approximation 

-BcsM = \lB1 -\- B'^ « B^{r csm), 

where Bt is a small transverse component. From several 3D 
mapping campaigns during which all coils (main and trim) were 
mapped, we know that this approximation is valid at the 10“"*^ 
level. 


4.1 Gradient extraction 


We extracted the vertical gradient by fitting a harmonic poly¬ 
nomial expansion of the magnetic field to the CsM array data. 
The choice of harmonic polynomials ensures that the result¬ 
ing expressions satisfy Maxwell’s equations. Due to the lim¬ 
ited number of magnetometers the expansion was limited to 
the next-to-linear order (NLO), which involves 9 parameters: 


Bz{x, y,z) =bo + gxX + gyy + gz z+ 

gxxix z ) -f gyy{y z )-t- 

gxyxy + gxzXz + gyzyz. (8) 

From expression ([^, one can easily calculate the volume av¬ 
erage of Bz and of its vertical gradient, assuming a trap with 
cylindrical symmetry: 


Bo = {Bz) 


bo + {gxx + gyy) 



dBz 

dz 


= 9z- 



(9) 

( 10 ) 


Let us now turn to the delicate task of estimating gradi¬ 
ent uncertainties. The two main sources that have to be taken 
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into account are the error on the magnetic field, and the ex¬ 
traction procedure. To assess their respective effects, extensive 
studies have been carried out using a toy model to generate 
known field distributions and check the extracted parameters 
f23i . It was found that the errors coming both from the mag¬ 
netometer offsets and from the expansion truncation never ex¬ 
ceed 5 pT/cm. In addition, we used a technique known as the 
jackknife method to get an error directly from the data. It in¬ 
volves performing a series of minimizations (unweighted in 
our case) by removing one out of the 16 magnetometers at a 
time. The dispersion of the extracted parameters provides an 
estimate of the error. For the different field configurations, we 
systematically obtained errors in the range 10 ± 5pT/cm, con¬ 
sistent with the model outcome. These jackknife errors were 
used subsequently in the analysis. 


4.2 Frequency shift measurement 

A sample of a raw data time series /ng against cycle number 
is displayed in Fig. together with the corresponding high- 
voltage values. Despite the large point-to-point fluctuations and 
a slow linear drift, one can clearly observe a small but system¬ 
atic correlation of the frequency shift with the electric field po¬ 
larity. To correct for the slow magnetic field drift, we sliced 
the data relative to the electric field polarity and analyzed data 
sets corresponding to the (-1- — —h) HV pattern. By doing 
so, any linear drift is exactly cancelled and higher orders are 
attenuated. 

For a data slice (-1-h), corresponding to 40 cycles, the 

extracted frequency shift and its uncertainty are given by 

^/ug = (/ng) ~ (■^Hg) 

and 

= y^(/4)' + ^(/ig)'. (12) 

where and stand for the mean frequency 

and its uncertainty as derived from the frequency distribution 
for the given HV polarity (-1- or —). Finally, a weighted mean 
over the whole set of data slices was performed to estimate 
the electric-field induced frequency shift 5fi^{E) for a given 
vertical gradient. 


5 Results and discussion 

The final result is displayed in Fig. The motional false mer¬ 
cury EDM is plotted against the extracted vertical gradient g^. 
The solid lines (red for sj, blue for Bq) correspond to a global 
linear fit with a single free parameter, namely the slope a ix^ jv = 
2.1/7). 

We can now compare the measured slope to its theoretical 
expectation from Eq. 6: 

loexpl = 1.122(35) X 10-2^ecm/(pT/cm), (13) 

andjothl = 1.148 x 10“^^ecm/(pT/cm). (14) 



Fig. 4. Mercury frequency versus cycle number. The blue line shows 
the value of the applied high voltage. 



Fig. 5. Motional false mercury EDM versus the vertical gradient 
for Sj (red up triangles) and S/ (blue down triangles). The solid lines 
correspond to a linear fit, and the dashed line to the theory discussed in 
section]^ The horizontal error bars are smaller than the symbol size. 


The agreement at the Icr level makes us confident that our 
magnetic gradient extraction procedure is reliable. This encour¬ 
aging result is nonetheless not sufficient to directly control the 
mercury-induced false neutron EDM at the required level of 
sensitivity. Indeed, for an error of 10 pT/cm on the vertical gra¬ 
dient, Eq. (|^ translates to a systematic error of 4.4 x 10“^®e cm 
on the neutron EDM, which is already larger than the current 
limit. There is fortunately a way to circumvent this issue. In 
their last nEDM paper fSl, the authors describe an analysis 
technique that enables one to find experimentally the work¬ 
ing point with no vertical gradient and therefore no motional 
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false EDM. This method, based on a tiny center of mass off¬ 
set between the cold neutrons and the warmer mercury atoms, 
nevertheless induces some additional systematic errors. These 
errors were carefully assessed and found subdominant with a 
final result statistically limited. 

Whereas the use of the Hg comagnetometer is essential 
and does not limit our nEDM sensitivity for the time being, 
with a foreseen sensitivity of a few 10“^^ecm in the coming 
years, new magnetometry solutions will be needed in the fu¬ 
ture. We pursue an intensive R&D program on magnetometry 
using ^^^Cs but also ^He atoms ll24l . In particular, efforts to¬ 
wards improving the absolute accuracy of the Cs magnetome¬ 
ters are currently underway ll25l as well as the implementation 
of Cs vector magnetometers ll26l . In parallel, we have started 
design and construction of a next generation nEDM spectrom¬ 
eter Q which, among other improvements, will benefit from 
a much better magnetic field control (passive and active). Ad¬ 
vanced magnetometry and improved magnetic shielding will 
be combined with co-magnetometry or could even allow opera¬ 
tion with only external magnetometers. Any possible mercury- 
induced false motional nEDM will therefore be much further 
suppressed or completely avoided. 


6 Conclusions 

We have performed a measurement of a frequency shift prmor- 
tional to the electric field strength for stored ^®®Hg atom^ us¬ 
ing a spectrometer devoted to the search for the neutron electric 
dipole moment at PSI. This shift, which we call the motional 
false EDM, originates from the combination of vertical mag¬ 
netic field gradients with the motional magnetic field and could 
be measured for the first time thanks to the unique combination 
of a mercury co-magnetometer and an array of external cesium 
magnetometers. The agreement with a prediction based on the 
general Redfield theory of relaxation provides additional con¬ 
fidence in the validity of our gradient-extraction procedure as 
well as in our capability to measure and control the vertical 
gradient. The same method was used in a recent measurement 
of the neutron magnetic moment lIZTll . 
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